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Electron scattering off an Anderson impurity immersed in the bulk of a 3D topological insulator is 
studied in the strong coupling regime, where the temperature T is lower than the Kondo temperature 
Tk- The system displays either a self-screened Kondo effect, or a Kondo effect with SO(3) or 
SO (4) dynamical symmetries. Low temperature Kondo scattering for systems with SO (3) symmetry 
displays the behavior of a singular Fermi liquid, an elusive property that so far has been observed 
only in tunneling experiments. This is demonstrated through the singular behavior as T —> 0 of 
the specific heat, magnetic susceptibility and impurity resistivity, that are calculated using well 
known (slightly adapted) conformal field theory techniques. Quite generally, the low temperature 
dependence of some of these observables displays a remarkable distinction between the SO(n=3,4) 
Kondo effect, compared with the standard SU(2) one. 

PACS numbers: 71.10.Pm, 73.43.-f, 72.15.Qm, 73.23.-b 


I. INTRODUCTION 

In this work we continue our study of a system com¬ 
posed of an Anderson impurity d immersed in a 3D topo¬ 
logical insulator (3DTI) with an “inverted-Mexican-hat” 
band dispersion around the F-poinlA— . For the sake of 
self-consistency let us very briefly recapitulate the pecu¬ 
liar feature of the pertinent system: Due to the special 
band structure, the ensuing Kondo effect (KE) is pro¬ 
foundly distinct from its metallic analog in normal met¬ 
als, because in addition to the original Anderson impurity 
d, there is an in-gap bound state (henceforth denoted as 
an / impurity), that is formed as a result of potential 
scatteringiii^. The d and / impurities form a “compos¬ 
ite quantum impurity” (CQI) that turns the pertinent 
Kondo physics much richer. The reason for that is as fol¬ 
lows: When isolated, the CQI can host two electrons that 
are found in singlet or triplet states with corresponding 
energies Es and Ep- Due to hybridization of the localized 
electrons in the CQI with the band electrons, the levels 
Es and Ep are renormalized with decreasing bandwidth 
albeit with different rate^. As a result, there is either 
a self-screened KE, or a KE with SO(3) or SO(4) dy¬ 
namical symmetry (depending on whether at the end of 
renormalization, Ep lies above, below or coincides with 
E;^. 

In our previous worb^, we have analyzed the pertinent 
Kondo physics in the weak coupling regime using pertur¬ 
bative RG analysis techniques. The goal of the present 
work is to perform an analysis of the Kondo scattering 
with SO(3) and SO(4) dynamical symmetry in the strong 
coupling regime T < Tk- The main physical motivation 
is to elucidate the occurrence of singular Fermi liquid in 
the SO(3) symmetric sector, that is exposed only in the 
strong eoupling regime. More than a decade ago it has 


been shown that the traditional classification of quantum 
impurity models into Fermi liquids and non Fermi liquids 
should be modified in such a way that one has to distin¬ 
guish between regular Fermi liquids and singular Fermi 
liquids (SFL)iiA^. The former case is exemplified by the 
standard SU(2) Kondo model where electrons are scat¬ 
tered from a magnetic impurity of spin <5=1 and at zero 
temperature the impurity is fully screened. This makes 
it possible to describe the system in terms of Nozieres 
Fermi-liquid picture. On the other hand, it was shown in 
Refii^ that when the impurity is under-screened, the cor¬ 
responding Fermi liquid is singular. Practically, it implies 
that the density of states is (logarithmically) singular at 
the Fermi energy. 

As is already stressed, this manifestation of SFL oc¬ 
curs only at low temperature T < Tk, and that re¬ 
quires the calculations to be carried out in the strong 
coupling regime, where perturbation theory is not ap¬ 
plicable. One need to resort to other approaches, such 
as Bethe ansatz or conformal held theory (CFT, that 
is employed here). Calculated experimental observables 
for SO(3) and SO(4) symmetric Kondo effect include the 
impurity contribution to the temperature dependence of 
the specific heat, magnetic susceptibility and resistivity. 
Comparing these results with those of the standard SU(2) 
Kondo effect, we indeed elucidate the remarkable SFL 
nature of the KE for the SO(3) symmetric sector. 

Experimentally, under-screened Kondo effect (USKE) 
has been observed in tunneling transport measurements 
through molecular transistorsiA2£. It is shown therein 
that the tunneling conductance G approaches the uni¬ 
tary limit when the temperature T tends to zero, but 
the derivative dG/dT is divergent. The possibility of de¬ 
tecting USKE in Kondo scattering in the bulk of metals 
is hampered by the fact that quantum impurities with 
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high spin give rise to Coqblin-Schrieffer type of Kondo 
scattering which eventually flows to a regular Fermi liq¬ 
uid fixed pointji^. In this work we overcome this obsta¬ 
cle and demonstrate the feasibility of observing USKE 
in bulk materials employing the occurrence of composite 
impurities. 

The organization of the paper is facilitated by the fact 
that the model and the relevant starting Kondo Hamil¬ 
tonian as well as the Kondo temperature have already 
been derived in our previous work (starting from the sin¬ 
gle impurity Anderson Hamiltonian). Therefore, in Sec. 
|H]we start right away by writing down the Kondo Hamil¬ 
tonian Hk describing low-energy exchange scattering of 
the band electrons by the d and / impurities as well as 
an exchange interactions between the d and / impurities. 
The Hamiltonian Hk possesses different dynamical sym¬ 
metries, SO(3) or SO(4), within various energy domains. 
The local density of states (DOS) of the TI with magnetic 
impurity immersed in it is calculated in Sec. IHII Kondo 
Interaction of the ” dressed” Fermi liquid (formed by the 
band electrons and the d impurity) with the /-impurity 
is discussed in Sec. IlYl In Sec. El we study and calculate 
the temperature dependence of the impurity contribution 
to the specihc heat, while the magnetic susceptibility of 
the impurity is considered in Sec. EH Sec. IVHI is devoted 
to the calculations of electric resistivity. The results are 
then summarized in Sec. Emi Some technical details 
are relegated to the Appendices. In Appendix lA II we 
describe the ground state of the isolated composite im¬ 
purity. Weak coupling renormalization of the exchange 
interaction strength of the /-impurity with the band is 
considered in Appendix lA 21 


II. KONDO HAMILTONIAN 

As it is shown in our previous worb^, the lowest-energy 
states of the isolated CQI are the singlet and triplet states 
with one electron on the d-level and a second electron on 
the /-level. Electron tunneling between the d-level and 
the band modifies the number of electrons in the CQI. 
Integrating out high energy states from the band edges 
renormalizes the singlet and triplet levels until charge 
fluctuations are quenched and one arrives at the local mo¬ 
ment regime. In this regime, the Schrieffer-Wolff trans¬ 
formation is used to map the Anderson Hamiltonian onto 
an effective Hamiltonian H = Hq -|- Hk- Here the first 
term, Hq, describes electrons in the bulk of the TI, 

Hq = (la) 

i/crk 

where v = ±1 denotes the conduction and valence band, 
efc = -I- {hvkY, Mk = mv^ — (lb) 


Mexican-hat” form with dispersion minimum at a surface 
of nonzero wave-vector q’s, with 



q = 


V 

hv 



( 2 ) 


where q = |q|. 

The Kondo Hamiltonian Hk assumes the following 
formal 


Hk = Jd{Sd ■ s) -t- Jf{Sf ■ s) -1- Jdf{Sd ■ S/), (3) 

where or Sj is the localized spin of the d- or /- 
impurity, s is the spin operator of the band electrons, 

L'kcr,i/''k'o’' 


T = (f®,f^,f^) is the vector of the Pauli matrices, 
couplings Jk, Jf and Jdf are explicitly given as^. 


Jd 

Jdf 


— id' 


PjAfil- 


Jf ^ 

‘JviPc 




f 




The 


(4a) 

(4b) 


where 


/3/ = 


V2 Vdf 
Af 


= ^f-^d + Uf. 


Here Cd and e/ are single electron energies of the d- and 
/-impurities, Vd is hybridization rate of the d-impurity 
and the band and V^f is the hybridization of the d- and 
/-impurities [see Fig. [5] below], and Uf is the Coulomb 
blockade parameter for the /-impurity. The Coulomb 
blockade Ud of the d-impurity is assume to be infinity^. 
Notice that, generally, d/ <C Jd^- 

The low-energy physics of the model depends on Jdf- 
More concretely, there are three different regimes for Jdf 
determining the different ground states (GS)^d£: 

• When Jdf = 0, the impurities are decoupled, 
and Kondo scattering of the band electrons is de¬ 
termined mostly by the hybridization of the d- 
impurity with the band, whereas the /-impurity 
can be considered as an isolated magnetic moment 
(i.e., it is not coupled to the band or to the d- 
impurity). In this case, there is standard Kondo ef¬ 
fect with full screening of the spin of the d-impurity. 
The exchange interaction of the /-impurity with 
the band manifests itself just at temperature below 
Kondo temperature. Therefore, Jf enters into con¬ 
sideration only for temperatures below the Kondo 
temperature. 


is the band dispersion. Accordingly, Sk is gapped and the 
insulator is topological for Bm > 0. For Bm >1/2 (as¬ 
sumed hereafter), the band dispersion has an “inverted- 


• When Jdf < 0 with \Jdf\ >> Jd, the ground state of 
the CQI is a triplet and the singlet state is highly 
excited. In this case, band electrons are scattered 
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from a CQI with spin S = 1, and there is an 
USKE of the magnetic impurity. Since the mag¬ 
netic moments of the d- and /-impurities in the 
triplet ground state are parallel to one another, the 
exchange interaction strength J/ of the /-impurity 
with the lead gives rise to slight modification of Jdr, 
(that is the exchange interaction strength of the d 
impurity with the band electrons). Therefore we 
assume J/ = 0 in this regime. 

• When Jdf > 0 with J^f ^ Jd, the ground state 
of the CQI is singlet and the triplet state has high 
excitation energy. In this case, there is no Kondo 
effect. 


e 



FIG. 1: (Color online) Dispersion e*,, Eq. Illbl l. ki {i = 1,2) are 
two solutions of the equation = tp- Ui = i^k- are the Fermi 
velocities at the inner or outer Fermi surfaces (fc = fci or ^ 2 ). 


Before presenting our calculations pertaining to physi¬ 
cal observables in our system, let us recapitulate the na¬ 
ture of the GS, specific heat and magnetic susceptibility 
for the ordinary case of a magnetic impurity immersed 
in an ordinary metal 0 . For the fully screened (FS) and 
the under screened (US) cases of the KF they are listed 
in the table below. Recall that the GS of the FSKF is a 
Fermi a liquid (FL) while that of the USKF turns out to 
be a SFL. Explicitly, 


the following expression for the phase shift, 


= 0 + 


1 , /£(^+i+lWM 

2* "l^r(^+i-^ln(^)) 


In 

2i 




^r(5+^lni 

For S = h, the above expression takes the form, 


( 6 ) 


KE 

GS 

Gy 

X 

FS 

FL 

- T 


US 

SFL 


^T-On-(^) 


In the following, we will consider physical properties of 
the system in the strong coupling regime (when the tem¬ 
perature is below the Kondo temperature) for the cases 
Jdf < 0 and Jdf = 0 , in turn. 


III. DENSITY OF STATES 


(5i(a;) = |-arctan(^), (7) 

whereas for 5* = 1 and uj <C Tk, <^ 5 ( 0 ;) has the asymptotic 
expression, 

Note that (5i(w) demonstrates a singular behavior near 
the point w = 0, whereas Si (uj) is regular. Differentiating 
([ 6 ]) gives the following expression for the DOS^^: 


In this section we apply the already known Bethe 
ansatz expression for the scattering phase shift and write 
down the density of states (DOS) that is peculiar to the 
system of spin S GQI immersed in a material with inverse 
Mexican hat band structure. Specifically, we focus on the 
low-energy spectrum of the system where it is justified to 
approximate the dispersion (llbll by linear expressions in 
k — ki {i = 1 , 2 ), where ki and ^2 are two solutions of 
the equation Sk = ep (see Fig. [J). The DOS Ns{oj) of 
the Fermi gas near the impurity position is expressed in 
terms of the phase shift Ss{uj) through the Friedel sum 
ruleH, 


^s(w) = -— Re 
ZtTU! 




(9) 


Here S' = 1/2 or 1 is the spin of the impurity, and the 
function /3(x) is defined as 



where is the digamma function, 

(ilnr(x) r'(x) 


?/'(x) = 


dx 


m' 


( 10 ) 


/ N 1 dSsiuj) ,,, 

Nsiuj) = - ( 5 ) 

TT diU 

where S = ^ or 1 is the impurity spin. Calculations 
based on the Bethe anzatsi^ applied for the full screened 
S = 5 and the under-screened S = 1 Kondo effect, yield 


The DOS Eq. (0 is shown in Fig. H for S = i and 
S = 1. It is seen that the DOS for S = 1 is singular. As 
a result the conventional Fermi liquid expansion of the 
phase shift can not be carried oulji^. The origin of this 
singularity is the non-analytic behavior of the phase shift 
Si{uj), eq. near the point w = 0 . 
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FIG. 2: (Color online) DOS Eq. lO for the SO(4) KE (blue curve) 
and the SO(3) KE (red curve). 


IV. KONDO INTERACTION OF THE FERMI 
LIQUID WITH THE /-IMPURITY 


When \Jdf\ Tk^, [where Tk^ is the Kondo temper¬ 
ature for the S0(4) KE, see Ref. — and Eq. (IA8I) in 
Appendix I A 2| . the d-impurity and the conduction band 
electrons form a singlet state, and the system composed 
of the d-impurity and the band electrons forms a local 
Fermi liquid that can be described within Nozieres Fermi 
liquid theoryii. The /-impurity is coupled to this local 
Fermi liquid through an effective exchange Hamiltonian, 

hP = Jf{Sfs), (11) 

where 

jf = MEisA, ( 12 ) 

Po 

The density of states of the local Fermi liquid is. 


P(e) = 




Ki 


(e - ep) + T|-^ 
The dimensionless coupling jf is. 


jf = JfPi^p) 


3j£q 

Tk, 


m ""p [T. 


(13) 


(14) 


where we take into account that pc ^ l/£g- The scaling 
equation for jf is. 


d\nD ^ 


(15) 


The initial value jfiTx^) of jf{D) is given by eq. (IT4l) . 
The solution of the scaling equation (US is, 


Jf{D) 


1 



(16) 


where the scaling invariant, the second Kondo tempera¬ 
ture, is 

exp I - — — exp —"j I. (17) 

y Jf \ JdJ j 

For Cd = —SOcq, Eq = 24eq, ep = 2eq, j = 0.08 and 
jf = 0.025, we get Tr^ = 0.15£g and = 0.002irR'^. 

Having set up the calculation framework we turn now 
to elucidate numerous physical observable. These include 

specific heat, magnetic susceptibility and impurity resis- 
( 2 ) 

tivity. Since is very small, we shall restrict our- 

( 2 ) 

selves to the temperature regime T ^ in the fol¬ 
lowing. For the SO(4) KE, the temperature is in the 
range Tr^ > T :s> T^\ For the SO(3) KE, we consider 
T < Tr^, where Tr^ and Tr^ are the Kondo tempera¬ 
tures for the SO(3) and SO(4) KE, respectively^. For the 
SO(3) KE, we employ CFT techniques, whereas for the 
SO(4) KE, we will apply CFT techniques for the inter¬ 
action between the d-impurity and the conduction band 
electrons, and then we employ poor man’s scaling for¬ 
malism to take into account the interaction between the 
local Fermi liquid and the /-impurity. 


V. SPECIFIC HEAT 


The first observable to be calculated is the specific heat 
of the TI within which the CQI is immersed. Its defini¬ 
tion reads. 


= ^/ deef{e)Ns{\e-ep\), 

where /(e) is the Fermi function. Noting that 

df{e) ^ e-ep 

dT ~ 4T2cosh2(^)’ 

we have 


a 


(0) 

V 



Ns{2Tx) dx 
cosh^(a;) 


(18) 


(19) 


For the SO (3) KE, the specific heat is Cy = Cy^ [given 
by Eq. m]- For the SO(4) KE, there is contribution to 
the specific heat due to the /-impurity. This contribution 
can be written as. 


37r2 

= — j1(T), (20) 

where jf{T) is given by Eq. (ITB)) . The total specific heat 

for the SO(4) KE is Cy = dy'’ + 5Cy. 

The ratio Cy jT [where Cy is the specific heat (fTOl) ] 
is shown in Fig. [3|for the SO(4) and SO(3) KE. In this 
case, the Kondo temperature is Tr^. For comparison, the 


In 
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FIG. 3: (Color online) Ratio Cy jT [where Cy is the specific heat 
1 II 9 I 1 I for the SO(4), SO(3) and SU(2) KE (solid blue, solid red and 
dashed blue curves). Here Tx = 7x4 SO(4) and SU(2) KE 

or Txg for the SO(3) KE. 


specific heat for the SU(2) KE ( i.e., for the case when 
the /-impurity is absent), is shown as well. It is seen that 
for the SU(2) KE, Cy jT saturates to constant as T —>• 0, 
in agreement with the prediction of Fermi liquid theory. 
The ratio Cy jT for the SO(4) KE is close to that for 
the SU(2) KE for T close to Tk^,- For low temperatures, 
Cy jT for the SO(4) and SU(2) KE deviate substantially 
from one another. This is result is due to the exchange 
interaction of the /-impurity with the local Fermi liquid. 
For the SO (3) KE, the ratio Cy jT diverges indicating 
the feature of a SEE behavior—. 


VI. MAGNETIC SUSCEPTIBILITY 


When Jfij < 0 and the ground state of the isolated CQI 
is S' = 1, the magnetic susceptibility is, 

x(T) = 2A.|ydeiVi(|e-e;^|) (21) 

where fVi(e) is given by Eq. ([5]) with S = 1. Taking into 
account that 


9/(e) ^ 1 

Se 4rcosh2(^)’ 


we can write, 


X{T) 


‘^XqTk /■ IVide-eFl) de 

T J cosh^ 


( 22 ) 


(23) 


where 


Xo = 



(24) 


For the case Jdf = 0 (see Eq. (g])), the susceptibility is 
calculated in the following way: First, we employ the 
GET technique to calculate the susceptibility for // = 0 
(see Eq. (flBll b Within this approximation, the CQI splits 
into two noninteracting impurities, one of them is coupled 
to the band electrons, and the other is fully decoupled. 


In this case, the susceptibility is given by. 


X^°\T) 


2nl J deiV.(le-eFl) 




(25) 


where Ni(e) is given by Eq. with S = ^. The sec¬ 
ond term, /r^/T, is the susceptibility of the /-impurity. 
Taking into account Eq. (0^ . we can write. 


x^°Ht) 


‘^XqTk 

T 


^i(k-eE|)de /i| 

cosh^ (^r 


(26) 


where xo is given by Eq. (1^ . The contribution to the 
susceptibility due to the interaction of the local Fermi 
liquid with the /-impurity is. 


Sx{T) = 


XqTk 

T 


hiT), 


(27) 


where jf(T) is given by Eq. (flBl) . The total susceptibility 
for the SO (4) KE is 

X(T) = x^°\T) + Sx{T). (28) 


The functions Tx(T) for the SO(3), SU(2) and SO(4) 


(TxWkXo) 



FIG. 4: (Color online) The quantity Tx(T) for the SO(3), SU(2) 
and the SO (4) KE is displayed as a function of temperature. The 
saturation as T ^ 0 for the SU(2) and SO (4) KE is reminiscent of 
the Curie law. On the other hand, for the SO(3) KE, Tx(T) 0. 


KE are shown in Fig. SI The SU(2) KE corresponds to 
the case when the /-impurity is absent. For this case, the 
Kondo temperature is Tk^- It is seen that for the SU(2) 
or SO(4) KE, Tx{T) TkXo for T ^ 0. For the SO(3) 
KE, Tx{T) vanishes as ln“^(T/Tx) when T —>■ 0, which 
is a manifestation of a singular Fermi liquid fixed point. 


VII. IMPURITY RESISTIVITY 

In this section we will calculate the impurity contribu¬ 
tion to the resistivity. Here, the peculiar band structure 
of the TI plays a central role. The calculation method 
of the impurity contribution to the resistivity depends 
on the underlying symmetry. The GET technique used 
above to derive the specific heat and magnetic field will 
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be employed to derive the resistivity of the KE with the 
SO(3) symmetry. On the other hand, for the KE with 
the SO(4) symmetry, we apply a somewhat different ap¬ 
proach: First, we apply the OFT technique to derive the 
resistivity for the case jf = 0 (see Eq. (fTHl) '). In the next 
step, we apply the perturbative RG (poor man’s scaling 
approach) to get the correction to the resistivity due to 
the interaction of the /-impurity with the local Fermi 
liquid. 

For jf = 0, the impurity resistivity can be written as^^ 

Here the conductivity (T) is 




3 


f d^k 

dfiek)' 

1 

dtk 


vlTtr{k), 


(29b) 


where /(e) is the Fermi-Dirac distribution, and Vk = |vk| 
is the group velocity, 

- i ^ - i k— 

h d\i. h dk 

For the dispersion (flbl) . Vk is given explicitly by, 

Vk = — {2Bm - l) , (30) 

£fe \q J mv 


where 


£o = mv^. 


The resistivity is directly related to the inverse of the 
transport relaxation time rtr(fc) that is expressible in 
terms of the phase shifts rji{k) (corresponding to angular 
moment 1) appearing in the partial wave expansion of the 
electron scattering from the impurity [see Eq. (2.20) in 
RefJ^]. Explicitly, 

1 

= (31) 


where 


Pc 


V 



ffi(e) 




i + (-i)' 



(33d) 

(33e) 


i = 1,2. Here 51 (e) and 52 (e) correspond to the two 
solutions of the equation efc = e. 

For pure s wave scattering, 50 (fc) = Ss{\ek — e^l) and 
5 / = 0 for / 7 ^ 0. Then Eq. (IHTl) can be written as. 


1 

Ttr(fc) 


2 Cinip 

Trhvik) 


sin^ 5s(|efe 


e^D- 


(34) 


The scattering phase i 55 (w) is given by Eq. ([ 6 |). For 
w = |e - cfI -C Tk, Ss{oj) is given by eq. © for 5 = i 
and by eq. ([ 8 |) for 5' = 1. 

Zero temperature resistivity: As T —0 we can write 


-^=S{e-ep). (35) 

Substituting Eq. (1551) into Eq. (I29b|) and taking into 
account Eqs. ® , © and dMl), we get 


r( 0 )( 0 ) = 


7re 


3n 


imp 




S{ek - ep) vl v{k). 


When the Fermi energy is constrained according to 


£q < ep < £ 0 , 


the equation £k = ep has two solutions, k = ki and k = 
^ 2 , (see Fig. [1}. The equation for the zero-temperature 
resistivity then takes the form. 


p(o)( 0 ) 

^^( 0 ) 


r(O)( 0 )’ 


a(o)( 0 )=ar( 0 ) + a: 


- 


(0) 




3ci, 


i = 1,2. 


( 0 ), 

(36) 


where p{k) is the bare density of states (DOS), 

Hk) = ek'), (32) 

k' 

(5(e) is the Dirac delta function and L is the linear size of 
the bulk. For the dispersion (llbl) . the DOS is 


v{k) = Fi(efc) -I- V2{ek), 

Pi(e) = i?(|e| - eg) d(eo - jej) 


Pc|e|5i(e) 



P 2 {e) = i?(|e| 


£ 9 ) 


Pc|e|52(e) 



(33a) 
, (33b) 

(33c) 


Finite temperature resistivity: When the tempera¬ 
ture is below the Kondo temperature Tk, the conductiv¬ 
ity (I29bl) can be written as. 


r( 0 ) 


iT) = x-y 


k^dk 


37r2 y 4Tcosh2(^i£=^’ 


vlnvik) + 


2T 


hrdk 


37r2 J 4Tcosh^{^i^) 


vlnr{k). 


Substitute into each of integrals a factor unity, 

J de S{e — £k)=l, and taking into account that for k < q 
or k > q, the equation £k = e has a single solution, fci(e) 
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or k 2 {e) (see Fig. [T]for illustration), we can write 


where jf{T) is given by Eq. (fT6|) . The constant Rq is, 


r( 0 ) 


(T) = 


3 


t^i(e)u^(e)rtr(fci(e))(ie 

4rcosh2(^) 


where 


3 


^2(e)u|(£)rtr(fc2(e))d£ 
4rcosh2(^ 


2T 


(37) 


k,{e) = qg^{e), * = 1 , 2 , 

and the functions gi{€) are given by Eq. (I33e() . Tak¬ 
ing into account Eqs. (IMl) and (I33ap . we can write the 
conductivity (l37t as. 


cr(o)(T) = 7 (£)**?(£) sin ^ {Ss{\e - eF\))de ^ 

4Tcosh2(l^) 

€.q 

^ TTfie^ r vl{e)vl{e) sin~^ (dg(|£ - £F|))d£ 

3cimp7 4Tcosh2(^) 

(38) 

The integrands on the right hand side of Eq. (1551) have 
a factor cosh~^( '^~j.^ ) which is equal to 1 for £ = cj’ and 
rapidly vanishes for e — ep ^ T. The behaviour of the 
other factors, I'fie) and vf{e) {i = 1 , 2 ), depends on the 
ratio {ep — eq)/T. When T <C ep — Eg, then vfie) and 
vf{e) change slowly within the interval \e — ep\ T and 
can be safely replaced by yf{ep) and vf{ep). In what 
follows, we will assume the inequality Tk <C £f ~ £q^ 
Then for T < Tp the conductivity (1551) takes the form. 


where 


Rq 


3*rCinip 1 
he^pl uj -I- ’ 


(42) 


= Vki, * = 1 , 2 , 

ki are two solutions of the equation Sk = ep, see Fig. [T] 
Hence, the resistivity for the SO(4) KE reads, 

p{T) = p(°HT)+6p{T), (43) 

where p^^\T) is given by Eq. (I40p . The resistivity (HOI) 


p(T)/p(0) 



FIG. 5: (Color online) Resistivity (I40II as a function of temperature 
for the SO(3) and SU(2) symmetries (solid red and dashed blue 
curves), and the resistivity 13511 for the SO(4) KE (solid blue curve). 

as a function of temperature is shown in Eigure [5] for 
the SO(4), SO(3) and SU(2) symmetries. Recall that as 
T <C Tp, the resistivity for the SU(2) symmetry is given 
by, 

/”>(!■) ,»d">( 0 )|l-^+o(^)|. (44) 




J dTcoshy^) 


(39) 


(The limits of the integration can safely be changed to 
=Fcx)). Einally, the impnrity resistivity is, 


p(°)(T) 


p(°)( 0 ) 


r sin {6s{\e\))de 

J 4Tcosh^(^) 


(40) 


This behaviour is typical for the Fermi liquid zero tem¬ 
perature fixed point. Deviation from Eq. (I44|) in the 
cnrves for the SU(2) and SO(4) KE is due to the Kondo 
interaction between the /-impurity and the local Fermi 
liquid. The temperature dependence of the resistivity 
for the SO(3) symmetry is rather distinct. Indeed, for 
T <C Tp, we can write 


p(o)(r) 


p(°)( 0 ) 



where p^°^(0) = 1/ (0) is the zero temperature impu¬ 

rity resistivity, see Eq. (1561) . 

Contribution of jf to the resistivity; The impurity 
resistivity for the KE with the SO(3) symmetry is given 
by Eq. (HiH) . That for the KE with SO(4) symmetry 
is contributed also by the interaction of the /-impurity 
with the local Eermi liquid. In order to derive this con¬ 
tribution, we apply the perturbation poor man’s scaling 
techniqne. 


+o(^ln-3(^^)^|. (45) 

The singular temperature dependence of the resistance 
is the result of the singular energy dependence of the 
scattering phase ©• 

VIII. CONCLUSIONS 


dp{T) 



The present work is motivated by the qnest to elu- 
(41) cidate the elusive SEE behavior within bulk materials. 
























Hamiltonian of the system is, 


This property, associated with USKE, has so far been 
observed only in electron tunneling experiments through 
quantum dots^^^— , where it is demonstrated that the 
tunneling conductance G approaches the unitary limit 
as r —>-0, but the derivative dG/dT diverges logarithmi¬ 
cally. Elucidating SEE in metals based on the USKE (e.g 
iron immersed in cupper) is not obvious because magnetic 
impurities with high spin give rise to Coqblin-Schrieffer 
type of scattering which falls on a regular Eermi liquid 
fixed poinlji^. Here we achieved this goal by analyzing 
Kondo scattering of electrons off an Anderson impurity 
in a 3D topological insulator with an “inverted Mexican 
hat” band dispersion. We used the fact that the inter¬ 
play between the Anderson impurity and its induced in¬ 
gap bound state results in self-screened Kondo effect or in 
the Kondo effect with SO(4) or SO(3) dynamical symme¬ 
tries. Using the conformal field theory technique, we have 
calculated the low temperature {T <C Tk) dependence 
of the specific heat, the magnetic susceptibility and the 
electric resistivity of the impurity for both screened and 
under-screened Kondo effect. Physical properties of the 
impurity for the under-screened case demonstrate zero 
temperature singularity which corresponds to a singu¬ 
lar Fermi liquid phase. In addition to the SEE behavior 
exposed here, this system exposes an interesting screen¬ 
ing mechanism where the impurity / is screened by the 
quasi-particles of the Fermi sea that is formed when the 
band electrons and the d impurity form a singlet state at 
T < Tk- 
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Appendix A: Anderson Model and Scaling Equations 

Part of the material presented in this appendix has al¬ 
ready been developed in our previous paper— where we 
analyzed the same system in the weak coupling regime. 
It is included here for the sake of self-consistency. In the 
first part we write down the bare Anderson model and 
specify the peculiarities resulting from the special form 
of the band structure. In the second part we elaborate 
on the scaling equations and RG flow for the SO (4) dy¬ 
namical symmetry. 


1. Anderson Model 

The system considered here is schematically displayed 
in Fig. [ 6 ^. It consists of a topological insulator whose en¬ 
ergy band has the “inverted mexican hat” structure, and 
an Anderson impurity immersed in its bulk (denoted as 
d-impurity). Potential scattering of electrons on the im¬ 
purity result in the formation a mid-gap localized energy 
level^, denoted hereafter as an /-impurity. The effective 


H = Ho + H, + Ht. (Al) 

Here the first term, Hq, is the Hamiltonian (ITa)) describ¬ 
ing electrons in the bulk of the TE 

The second term on the right hand side of eq. (EH) is 
the Hamiltonian of the isolated CQI, composed of the d- 
and /-levels, 

H, = Hd + Hf + Hdf. (A2a) 

Here Hd and Hf are the (atomic) Hamiltonians of the 
d- and /-impurities, and Hdf describes the hybridization 
between them. 


Hd 

— ^ ^ '^dcr “ 1 “ 

a 

(A2b) 

Hf 

(j 

(A2c) 

Hdf 

= Ud/^(/K + 4/.), 

(A2d) 


<7 


where or e/ is the d- or /-impurity energy level and 
Ud or Uf is the interaction between electrons on the im¬ 
purity. Here Udo- = d\dcr, Ufa = flfa, d\ or da is the 
creation or annihilation operator of electron on the d- 
level, /j! or fa is the creation or annihilation operator of 
electron on the /-level. 

The last term on the right hand side of eq. lAH), Ht, is 
the hybridization between the Anderson impurity d and 
the band electrons, 

E + 47.k.) . (A3) 

k,1^,(7 

Note that hybridization prevails only between the band 
electrons and the d impurity. The /-level is formed as a 
result of the potential scatering of the conduction band 
electron on the d-impurity^. Therefore, hybridization of 
/-impurity and the band electrons is absent. 

Energy scales: Few words about energy scales are in 
order: Unless otherwise specified, we shall assume that 
Ud ^ oo and 

ef — Do < Cd <C e/ < ef < e/ + U/ <C ef + Dq, (A4) 

where Dq (the initial bandwidth) is the highest energy 
cutoff, and ef is the Fermi energy (see Figure [6]). In the 
calculations below we use ef = 2 eq, Cf Uf = bSq 

and Ed = —SOCg. 

The eigenstates of He, equation (jA2all . are specified 
by the configuration numbers {Nd,Nf) indicating the 
number of electrons on the levels d and /. With en¬ 
ergy scales specified in Eq. (Kill . the ground state (GS) 
has Nd = Nf = 1 and there are four possible states, 
a spin-singlet state \S) and three spin-triplet states \Tm) 
(to = 0, ±1). The singlet energy is modified when Vdf ^ 0 
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FIG. 6: (Color online) Panel (a): Illustration of the Anderson 
model consisting of a band, d- and /-impurities. The hybridization 
rates of the d-impurity with the band and the /-impurity are and 
respectively. Panel (b) Energy dispersion (|lbl> . and energy 
levels Ed and e f . Here e p denotes the Fermi energy. All the energy 
levels of the band below the Fermi energy are occupied (green area), 
whereas the white area above ep denotes unoccupied levels. 


while the triplet energy is unaffected. Explicitly, 

1^) = {^{4ft - 44) - A5a) 

|ri) = 4/t|0), |r_i) = d|/||0), 

(A5b) 

\To) = ^{4fl+4fl}\0), 


and jf = JfPo must proceed to third order in these 
parametersi^, that is, 


djd 

dlnD 

9jf 

dlnD 


—Jd + jd{jd + J/)> 
-j}+jf{jd +jf)- 


(A6a) 

(A6b) 


The initial values of jd and jf at D — Da are JdPc and 
Jf Pc, where Jd and Jf are given by Eq. (Hal) . Scaling of 


A 



FIG. 7: (Color online) Scaling of and jf. 

jd and jf is shown in the flow diagram. Fig. [71 The red 
dot denotes the two-channel fixed point j^ = jf = 1/2. 
We are interested in the solution of these equations under 
the following inequalities: 

jfi^ii) ^ jdiDii) <C 1, 


where 


£S = f-d + e/ 



as = ^1-/3/, Pf 

- ^d + Uf. 


et = Cd + e/, 


V2Vdf 


In the absence of hybridization of d-electron with the 
band electrons (i.e., when Vd = 0), eg < ex'- As expected, 
the singlet state has lower energy than the triplet state. 


2. Scaling Equations for the SO(4) Dynamical 
Symmetry 

In order to perform the RG analysis considered in 
Sec. HV] [see Eq. mi we need the effective coupling 
jf{TKi) at the effective bandwidth equal to Tk^- To de¬ 
rive jf {Txi ), we apply the weak coupling RG analysis for 
the SO(4) KE. In the the weak coupling regime, the scal¬ 
ing equations for the dimensionless couplings jd = JdPo 


so that we are far away from the two-channel fixed point. 
In this case jd{D) increases when D decreases, whereas 
jf{D) decreases and goes to zero when D vanishes. The 
solution of the scaling equation (IA6al) can be approxi¬ 
mated as. 


Jd[D) 



(A7) 


where the scaling invariant, the Kondo temperature Tx 4 , 
is 

7^4 = D^i exp I - I. (A8) 

The solution (IAtI) is valid just for D Tk^- A more 
satisfactory approximation is given in RefJ^, 


Jd{D) 



1 - 




In 




Stt^ 


(A9) 
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The solution of the scaling equation (IA6bl) is, 

^ Du ^ 


3 AD) 

jf{Dii) 


= expj -y^ * (AlO) 

^ D ^ 


Taking into account eq. dSH), we get 


JfjP) 

3f{Dii) 


= exp <^ — ( - T{D) 


where 


D 


HD) = JlA — + —-In — 


Tk. 


Stt^ 


D 


Tk, 


When D approaches Tk,, jfiTx,) can be approximated 
as 


jf{TK,) = jf{Du) exp II 

~ 0-Af{Dzi), 


(All) 


where we take into account that Da ^ Tk, ■ 
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